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Abstract. Wc study equal weight numerical integration, or Quasi Monte Carlo (QMC) 
rules, for functions in a Sobolev space H s (§ c! ) with smoothness parameter s > d/2 defined 
over the unit sphere E> d in R d+1 . Focusing on A-point configurations that achieve optimal 
order QMC error bounds (as is the case for efficient spherical designs), we are led to 
introduce the concept of QMC designs: these are sequences of A-point configurations Xm 
on E> d such that the worst-case error satisfies 



=5 



sup 

/eH s (s d ), 
II/IIh»<i 



0(N- s/d ), A 



00, 



with an implied constant that depends on the H s (S d )-norm, but is independent of A. 
Here aa is the normalized surface measure on S d . 

We provide methods for generation and numerical testing of QMC designs. An essen- 
tial tool is an expression for the worst-case error in terms of a reproducing kernel for the 
space H s (S d ) with s > d/2. As a consequence of this and a recent result of Bondarcnko 
et al. on the existence of spherical designs with appropriate number of points, we show 
that minimizers of the A-point energy for this kernel form a sequence of QMC designs 
for HF(§ d ). Furthermore, without appealing to the Bondarcnko et al. result, we prove 
that point sets that maximize the sum of suitable powers of the Euclidean distance be- 
tween pairs of points form a sequence of QMC designs for H s (S d ) with s in the interval 
(d/2, d/2 + 1). For such spaces there exist reproducing kernels with simple closed forms 
that are useful for numerical testing of optimal order Quasi Monte Carlo integration. 

Numerical experiments suggest that many familiar sequences of point sets on the sphere 
(equal area points, spiral points, minimal [Coulomb or logarithmic] energy points, and 
Fekete points) are QMC designs for appropriate values of s. For comparison purposes 
we show that configurations of random points that are independently and uniformly dis- 
tributed on the sphere do not constitute QMC designs for any s > d/2. 

If (Xn) is a sequence of QMC designs for HP(§ d ), we prove that it is also a sequence of 
QMC designs for M s (E> d ) for all s' € (d/2, s). This leads to the question of determining 
the suprcmum of such s, for which we provide estimates based on computations for the 
aforementioned sequences. 
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1. Introduction 

In this paper we introduce a new notion for sequences of finite point sets on the unit 
sphere § d in the Euclidean space IR d+1 , d > 2, namely that of sequences of QMC designs. 
These are sequences that emulate spherical designs in that they provide optimal order 
equal weight numerical integration (or Quasi Monte Carlo) rules for certain Sobolev spaces 
of functions over the unit sphere E> d . 

A spherical t-design, a concept introduced in the groundbreaking paper [22] by Delsarte, 
Goethals and Seidel, is a finite subset Xn C S d with the characterizing property that an 
equal weight integration rule with nodes from Xn integrates exactly all polynomials P 
with degree < t; that is, 



N 



degP < t. 



Here N = \X^\ is the cardinality of X^, or the number of points of the spherical design, 
while the integral is with respect to the normalized surface measure on E> d , and the 
polynomials of degree < t are the restrictions to S d of the polynomials of degree < t on 
R d+1 . 

Sequences of spherical designs have a known fast- convergence property in Sobolev spaces. 
(See Section l2~31 below for the definition of the Sobolev space IF(S d ).) This property, stated 
in the following theorem, was first proved for the particular case s = 3/2 and d = 2 in |27j . 
then extended to all s > 1 for d = 2 in [28] , and finally extended to all s > d/2 and all d > 2 
in [IH]. (The results in those papers were proved for all positive- weight integration rules 
with an appropriate degree of polynomial accuracy in relation to the number of points, but 
here we restrict our attention to equal weight rules.) 

Theorem 1. Given s > d/2, there exists C(s,d) > depending on the M s (E> d )-norm such 
that for every N -point spherical t-design X^ on E> d there holds 



(2) sup 

f£W B (§ d ), 

ll/llw«<i 



N ^ N hd 



< C(s,d) 



Note that the constant C(s, d) in this theorem does not depend on t or on N, or on the 
particular spherical design Xn- Note too that the condition s > d/2 is a natural one, since 
by the Sobolev embedding theorem this is the condition needed for H s (§ a! ) to be embedded 
in C(S d ). 

The relation between N and t in a spherical design is not fixed, but there are known 
lower bounds on N (see (jSJ) below) that tell us that N is at least of order t d , and a recent 
result [TTJ of Bondarenko et al. (see Theorem [H] below) asserts that given t there always 
exists a spherical design with N x t d . Here we write a n x b n to mean that there exist 
positive constants c\ and C2 independent of n such that c\a n < b n < C2a n for all n. 

Motivated by these facts and the belief that the only interesting sequences of spherical 
designs are those with N x t d , we now define the notion of a sequence of QMC designs. 
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Definition 2. Given s > d/2, a sequence (Xn) of Appoint configurations on S d with 
N — > oo is said to be a sequence of QMC designs for M. s (S d ) if there exists c(s,d) > 0, 
independent of N, such that 



(3) sup 

/GH s (S d ), 

II/IIef<i 



c(s,rf) 



In this definition Xn need not be defined for all natural numbers N: it is sufficient that 
Xn exists for an infinite subset of the natural numbers. By a special case of theorems in 
[26] and [25], the exponent of N in ([3]) cannot be larger than s/d: 

Theorem 3. Given s > d/2, there exists d(s, d) > depending on the M. s (E> d )-norm such 
that for any N -point configuration on E> d 



( 4 ) ~^r< sup 

/GH s (S d ), 



c'(s, d) 
N s / d 

j 

I/IIh«<i 



The following theorem, obtained by appealing to results of Brandolini et al. [13], asserts 
that if (Xn) is a sequence of QMC designs for H s (S d ), then it is also so for all coarser 
Sobolev spaces U s '(§ d ) with d/2 < s' < s. Like all statements in the paper needing formal 
proof, it is established in Section |9j 

Theorem 4. Given s > d/2, let (Xjy) be a sequence of QMC designs for HP(S d ). Then 
(Xn) is a sequence of QMC designs for IF (S d ) , for all s' satisfying d/2 < s' < s. 

It follows from this theorem that for every sequence of QMC designs (Xn) there is some 
number s* such that (Xn) is a sequence of QMC designs for all s satisfying d/2 < s < s* , 
and is not a QMC design for s > s*; that is 

(5) s*:=s*[(Xjv)]'.= sup {s : (Xn) is a sequence of QMC designs for H s (S d )}. 
If s* = +oo, we append the adjective "generic." 

Definition 5. A sequence of TV-point configurations (X N ) on B> d is said to be a sequence 
of generic QMC designs if (|3]) holds for all s > d/2. 

As in Definition [21 Xn need not be defined for all natural numbers N. Obviously, every 
sequence of spherical t-designs with x t d as t — > oo is a sequence of generic QMC 
designs for HP(S d ), for all s > d/2. We noted already the claimed existence of a sequence 
of spherical t-designs with x t d . A simple application of that result yields the following. 

Theorem 6. There exist N -point spherical t-designs Yn on E> d for N = 1,2,3,... and 
t x A^ 1 ^ that form a sequence of generic QMC designs. 

For fixed s > d/2, there exist many sequences of QMC designs for HP(S d ) that are 
not composed of spherical designs. Indeed, if K is a reproducing kernel for the Sobolev 
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space H s (S d ), s > d/2, we prove in Section [3] (by appealing to Theorem [6]) that iV-point 
configurations (N > 2) minimizing the i^-energy functional 



N N 



(6) 



j=i i=i 



form a sequence of QMC designs for this HP(§ d ); cf. Theorem [T2j For s in the interval 
(d/2, d/2 + 1) we show in Section [5] that, for C a suitably large constant, C — |x — y\ 2s ~ d 
is a reproducing kernel for IF(8 d ), and therefore the maximizers of the generalized sum of 
distances 



(7) 



N N 

EE I*, 

j=i t=i 



1 2s— d 



iV = 2, 3, 4, . . . form a sequence of QMC designs for this IfP(§ d ). 

We also provide an alternative sufficient condition for QMC designs that utilizes poly- 
nomial truncations of a zonal reproducing kernel K(pc, y) = K(x ■ y), but requires also a 
regularity condition (Property R) imposed on the point configurations; see Definition [14" 
and Theorem [T5l 

Numerical evidence presented later in this paper suggests that many familiar sequences 
of point sets on § 2 (such as minimal [Coulomb or logarithmic] energy points, generalized 
spiral points, equal area points, and Fekete points) form sequences of QMC designs for 
IfP(§ 2 ) for values of s up to a supremum s* that depends on the particular sequence. Some 
conjectured values of s* are given in Section El 

That the QMC design property ([3]) is not satisfied by all sequences of point sets follows 
from a probabilistic argument. 

Theorem 7. Given s > d/2, the expected value of the squared worst-case error satisfies 



\ 



E 



sup 

/GH s (S d ), 

II/IIhs<i 



i N r 

*£/w-/, /Wd 



2-i 



b(s, d) 



for some constant b(s, d) > 0, where the points Xi, . . . , Xjy are independently and uniformly 
distributed on E> d . 

Theorem [7] tells us that randomly chosen point sets give a slower rate of convergence 
than N~ s / d for all s > d/2, and hence do not form QMC designs. (See Section [7] for a 
more complete discussion.) However, if we compartmentalize the random point selection 
process with respect to a partition of the sphere into iV equal area regions with small 
diameter, then we do get an average worst-case error rate appropriate to QMC designs for 
s G (d/2, d/2 + 1), see Theorem 1241 On the other hand, such randomized equal area point 
configurations will, on average, not form a sequence of QMC designs for s > d/2 + 1, see 
Theorem |25j 

We shall also discuss "low-discrepancy sequences" on the sphere and estimates for their 
worst-case error when used for QMC rules. It turns out that the point sets of such a 
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sequence almost satisfy the QMC design property for s G (d/2,(d+ l)/2), except for a 
power of logiV. 

The structure of the paper is as follows: The next section provides background for 
spherical designs and for Sobolev spaces H s (S d ) and their associated reproducing kernels. 
Section [3] characterizes the worst-case error for equal weight numerical integration and 
states two main results: Theorems [12] and [15j Examples of Sobolev spaces and associated 
kernels are given in SectionsH]and[5j with particular emphasis on configurations maximizing 
sums of generalized distances. Section [6] concerns low-discrepancy sequences on the sphere 
and their quadrature properties. In Section [7] we analyze the quadrature error for randomly 
chosen points on the sphere and, in Section [SJ we provide numerical results for worst-case 
errors and quadrature errors for certain familiar sequences of configurations. Most of the 
formal proofs are given in Section |HJ 



2. Background 

2.1. Spherical designs. In the literature on spherical designs, and again in this paper, 
the relation between N and t in ([1]) plays an important role. It is known (Seymour and 
Zaslavsky [38]), that a spherical t-design always exists if N is sufficiently large but that 
result says nothing about the size of N. In the important paper [22] lower bounds of exact 
order t d were established, precisely 



N > 



'd + t/2\ fd + t/2-l\ 

d +[ d ] forteven, 



'd+\t/2\\ 
d J 



for t odd, 



but it is known, see Bannai and Damerell [51 Ej, that these lower bounds can be achieved 
only for a few small values of t. Korevaar and Meyers [3U] conjectured that there always 
exist spherical t-designs with N x t d points. Bondarenko, Radchenko and Viazovska [IT] 
claim to have resolved this long-standing open problem, by establishing the following result. 

Theorem 8. Ford > 2, there exists a constant Cd depending only on d such that for every 
N > c d t d there exists a spherical t-design on E> d with N points. 



Supporting evidence is provided in [18], which used interval analysis to establish rigor- 
ously the existence of spherical designs on § 2 with N = (t + l) 2 for all values of t < 100. 

2.2. Spherical harmonics. Recall that a d denotes the normalized (Hausdorff) surface 
measure on the unit sphere S d in R d+1 . The [ non-normalized] surface area of E> d is denoted 
by oJd- For future reference, we record the following facts: 

rn , iwd-i uj d ^ r((d+i)/2) d 1 ' 2 

(9) 7d-=-i , = ' ; ' : ~ — = as d -> oo. 
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Here, denotes the gamma function and f(x) ~ g(x) as x — > c means f(x)/g(x) — > 1 
as x — > c. The asymptotic relation in OH]) follows from (see [U Eq. 5.11.12]) 

(10) — — ~ z a ~ h as z — > oo in the sector I argzl < ix — 5 < it. 

T(z + b) 

We make use of the rising factorial, that is the Pochhammer symbol defined by 

(11) (a) = l, (a) n+1 = {a) n {n + a), n = 0,l,..., 

which can be written in terms of the gamma function by means of (a) n = T(n + a)/ 17(a). 

We denote, as usual, by {Y^ : k = 1, . . . , Z(d,£)} a collection of L 2 -orthonormal real 
spherical harmonics (homogeneous harmonic polynomials in d + 1 variables restricted to 
E> d ) of exact degree £, where 

(12) Z(d, n) = (2» + d - 1) ~ ^ Mn ^oo. 

It is well-known that the satisfy the following identity known as the addition theorem: 

Z(d,£) 

(13) £ y /ifc (x)y /jfc (y) = Z(d,£)P^\^- y), x,y e § d , 
fc=i 

where is the normalized Gegenbauer (or Legendre) polynomial, orthogonal on the inter- 
val [—1, 1] with respect to the weight function (1 — t 2 ) d / 2_1 , and normalized by P\ > (1) = 1. 
Each spherical harmonic Y^. of exact degree i is an eigenfunction of the negative Laplace- 
Beltrami operator — A d for S d with eigenvalue 

(14) Xf.=£ {£ + d - 1) . 

(For further details see, e.g., [55].) 

The family {Yp^ '■ k = 1, . . . , Z(d, £); £ = 0,1, . . .} forms a complete orthonormal (with 
respect to ad) system for the Hilbert space L 2 (S d ) of square- integrable functions on S d 
endowed with the usual inner product and induced norm 



(f,g)h 2 (&y-= / /(x)^(x)da d (x), \\f\\u{&)'=\/Uj) 



L 2 (S d )- 



We shall denote by P t (§ d ) the space of all spherical polynomials of degree < t (that is, 
the restriction to § d of all polynomials in d + 1 real variables of degree < t). The space 
P t (S d ) coincides with the span of all spherical harmonics up to (and including) degree t, 
and its dimension is Z(d + l,t). 

We make frequent use of the following simple application of the addition theorem. 

Lemma 9. Let d > 2. For all integers £ > and all choices of points xi, . . . , x^ G § d 
there holds 

N N Z{d,£) N 2 

(15) * i (x u ...,x„):=—^^z(d,£)pM(x j -x i ) = E nH y ^) 

j=l i=l k=l j=l 
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2.3. Sobolev spaces. The Sobolev space H s (S d ) may be defined for s > as the set of 
all functions / £ L 2 (S d ) whose Laplace- Fourier coefficients 



(16) 
satisfy 

(17) 



fe,k '■— (/> Y t 



oo Z{d,l) 
1=0 k=l 



/(x)Y^(x)da d (x) 



f* 



< oo, 



where the A/s are given in Q14p . On setting s = we recover H°(S d ) = L2(S d ). 

The norm in IF(§ d ) may of course be defined as the square root of the expression on 
the left-hand side of the last inequality. In this paper we shall, however, take advantage of 
the freedom to define equivalent Sobolev space norms. Let s > d/2 be fixed and suppose 

we are given a sequence of positive real numbers (o^)i>o satisfying 

(is) 4 s) x(i + A f r s x(i+ir 2s . 

Then we can define a norm in H s (S d ) by 



(19) 



oo z(d,e) 

E E ~Js) 

e=o k=i a e 



-| 1/2 



fu 



The norm therefore depends on the particular choice of the sequence (ty )e>o, but for 
notational simplicity we shall generally not show this dependence explicitly. Clearly, Def- 
initions [2] and [5] are not tied to a particular Sobolev norm || • since a change to an 
equivalent Sobolev norm merely leads to a change of the constant c(s, d). The correspond- 
ing inner product in the Sobolev space is 

oo z(d,e) 

(20) (/,<?V:=E E 



\k ge,k- 



=0 k=l 



It is well known that HP(§ d ) C H s '(§ d ) whenever s > s', and that H s (S d ) is embedded in 
the space of fc-times continuously differentiable functions C fc (S d ) if s > k + d/2 (e.g. |25j). 

2.4. Sobolev spaces as reproducing kernel Hilbert spaces. Since the point-evaluation 
functional is bounded in HP(§ d ) whenever s > d/2, the Riesz representation theorem as- 
sures the existence of a reproducing kernel K^ s \x, y), which can be written as 

oo oo Z(d,£) 

(21) # (S W) = ^4 s) Z(d,£)Pf (x-y) = E ^YeMYM, 

1=0 1=0 k=l 

where the positive coefficients a\ satisfy ( II8|) . It is easily verified that the above expression 
has the reproducing kernel properties 

(22) K {s \-^) EM s (S d ), x£§ d , (/,^(-,x)) Hs = /(x), x£§ d ,/£ff(§ d ). 
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The kernel is a zonal function; that is (x, y ) depends only on the inner product x ■ y . 
We write for simplicity K^(x ■ y):=A^(x, y). For the particular choice of = (1 + Xe)~ s , 
we use the notation 

oo 

(23) JifW (x, y) := + \t)- s Z(d, £)P«(x • y), 

1=0 

which we call the canonical kernel for H s (S d ). 

Sections S] and [5] contain explicit examples of Sobolev spaces and associated kernels. 

3. Numerical integration and worst-case error 

3.1. Worst-case error. Our results are based on an explicit expression for the "worst-case 
error" that occurs on the left-hand side of 

Definition 10. For a Banach space B of continuous functions on § d with norm || ■ ||#, 
the worst-case error for the integration rule Qf^yv] with node set = {xi, . . . , x^} 
approximating the integral 1(f), with Q[Xjv](/) and /(/) defined by 

i - r 

(24) Q[X N ](f):=-J2f(^), 1(f)- /(x)d^(x), 
is given by 

(25) wce(Q[X N ]-,B):=sup{\Q[X N ](f)-I(f)\ : / £ B, \\f\\ B < l} . 

As a trivial consequence of the definition we have the following error bound for an 
arbitrary function / 6 B: 

(26) \Q[X N ](f) - I(f)\ < wce(Q[X N ]; B) \\f\\ B . 

Because of the similarity of f|2^|) to the celebrated Koksma-Hlawka inequality, which in- 
volves the "star-discrepancy" of the node set, the worst-case error is sometimes referred to 
as a generalized discrepancy, see for example [21]. In this paper, however, we shall generally 
reserve the word "discrepancy" for quantities that have a geometric interpretation. 

3.2. Worst-case error in a reproducing kernel Hilbert space. For most reproducing 
kernel Hilbert spaces there is a computable expression for the worst-case error, as shown 
by the following standard argument. With K the kernel of a reproducing kernel Hilbert 
space H with inner product (■, •)#, the reproducing kernel property /(x) = (f,K(-,x)) H 
allows us to write 

Q[X N ](f) - 1(f) = (f,K[X N ]) H , feH, 
where 7^[Xat] £ H is the "representer" of the error, given by 

1 N 

n[x N ](xy.=- k(x, x,) - / y A'(x, ■), 

i=i 
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assuming that the integration functional / t— > 1(f) is bounded on H. Here I y K means 
the integral functional / applied to the second variable in K (and later I X K will mean the 
integral functional applied to the first variable) E 
It follows that 

[wce(Q[X N ];H)] 2 :=[sn P {\Q[X N ](f) - I(f)\ : f E H, \\f\\ H < l}'* 

sup{\(f,K[X N }) H \ :f e H, \\f\\ H < 1} 

= \\n[x N }\\ 2 H = (n[x N },n[x N }) H 

TV TV TV 
j=l j=l j=l 

3.3. Worst-case error in H s (§ a! ). Now consider the special case of the reproducing kernel 
Hilbert space H s (S d ) with s > d/2, and with reproducing kernel given by ( 121]) . For this 
case it is easily seen that 



from which it follows that 



(27) [wce(Q[X w ];H s (S d ))] ; 



/ y i^)(x, 



TV TV 



1 {S) 



j=l i=l 



TV TV 



-a, 



j=i t=i 



where /C (s) : [—1,1] ->■ K is defined by 

oo 

(28) /C^(x ■ y):= £ af Z(d, *)P«(x ■ y). 



£=1 



The use of the calligraphic symbol here and for subsequent kernels indicates that the sum 
runs from £ = 1 rather than I = 0. Note that the kernels depend on s through the sequence 

(a^ s) )£>o- 

We summarize these observations in the following proposition. 

Proposition 11. For s > d/2, let IFF (3 d ) be the Hilbert space with norm (1191) . where the 
sequence (of )e>o satisfies ffTSl) . and let K,^ be given by fT2"8l . Then, for a rule Q[X N ] with 
node set X^ = {xi, . . . , xtv} C E> d , 

1/2 



(29) 



/ TV TV \ 

wce(Q[X N }; W (§ d )) = — "£Y1 " x *) 

V i=i i=i J 



Z(d,£) 



EE 

=1 k=l 



1 * 
7^E y < 



fcl x T 



3=1 



1/2 



*The norm ||72.[Xjv]||if is also known as the g-diaphony of Xn with 5 = I y K(yi, •); see [3]. 
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From the first expression in ([29]) . the squared worst-case error for the rule QfXjy] is the 
normalized /C^-energy functional evaluated at the node set X^. This expression can be 
computationally useful when the kernel KS a > is available in closed form. 

By comparison with a sequence satisfying Theorem [6] we deduce that the minimizers of 
([29]) yield a sequence of QMC designs for H s (§ d ). 

Theorem 12. Under the assumptions of Proposition \TT\ if X* N , N — 2,3,4, ... , minimizes 
the energy functional 

N N 

j=i t=i 

then there exists c(s, d) > depending on the M. s (S d )-norm such that for all N > 2 

wce(Qra;ff(S d ))<^. 

Consequently, (X^) is a sequence of QMC designs for HP . 

The next result says, in essence, that the computed worst-case errors in H s (S q! ) of a 
given sequence of QMC designs for HP where d/2 < s' < s, will always show a rate of 
decay at least 0{N~ s '/ d ) and at most 0{N" s / d ). 

Theorem 13. Under the assumptions of Proposition [771 for HPfS ') and HP' (E> d ) , if (X^) 
is a sequence of QMC designs for HP (S d ), then 

(30) wce(Q[X N }; W (§ d )) < ^ffi , d/2 < s' < s, 

where c(s, s', d) > depends on the norms for HP (S d ) and H s '(S d ) 7 6itt zs independent of N. 

This result follows from the last expression in ( 129]) . since it implies that, with respect to 
the canonical kernels ( 123|) for HP(§ d ) and HP there holds for any iV-point configuration 

(31) wce(Q[X N }; ff(§ rf )) < wce{Q[X N ]; H s '(§ d )), d/2 < s' < s. 

We shall exploit Theorem [TBI in Section [S] to determine empirical values of the supremum 
s* in ([5]) for a number of putative sequences of QMC designs. 

3.4. Truncations of the Laplace-Fourier series. Using truncations of the Laplace- 
Fourier series of a fixed reproducing kernel, we provide in the next theorem a sufficient 
condition for a sequence of point sets to be a QMC design sequence. It relies on the 
following regularity condition, first introduced by Sloan and Womersley in [39] in the 
context of positive weight quadrature rules. This imposed condition allows us to control 
the contribution to the worst-case error of the suppressed tail-part of the aforementioned 
kernel. 
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Definition 14 (Property R). A sequence (Z^) of iV-point configurations on E> d is said to 
have the Property R if there exist positive numbers Cq and c\ independent of N, such that 
for all Zn and all z £ E> d the nodes satisfy 

\Z N r\C{z- Cl N~ l l d )\ <c , 

where 

(32) C(z;6):={y £ § d : y z > cosO} 

denotes the spherical cap of geodesic radius 9 centered at z. 

The regularity Property R expresses a natural requirement, that the number of points in 
a spherical cap with radius of order A _1 / d , and hence whose area is of order A -1 , should 
be bounded independently of both A and the location of the cap. 

Remark. A sequence of Appoint sets (Xn), Xn = {x^jv, • • • , X-n,n}, is called well- separated 
if there exists a constant q > such that 

mm |xjyv - Xi tN \ > N = 2, 3, 4, ... . 

It is easily seen that a well-separated sequence of configurations satisfies Property R, but 
not conversely. Furthermore, Reimer [35] has shown that every sequence of spherical t- 
designs with N x t d points automatically satisfies Property R. 

Theorem 15. Under the assumptions of Proposition [77], a sequence (X^) of N -point 
configurations on E> d satisfying Property R is a sequence of QMC designs for W(S d ) if and 
only if for integers t with t x N 1 ^ 

1 N N 

( 33 ) Jfi E E ■ = 0{N~^ d ) as N ^ oo, 

j=i i=i 

where K,\ is the truncated kernel corresponding to ([28]) . 

t 

(34) /Cj s) (x • y):= 4 S) Z(d, t)Pi d \x ■ y). 

l=\ 

This fits well with the framework of variational characterization of spherical designs (see 
[T9] [20] |24"[ |4T|). where one seeks point configurations for which the left-hand side in ([33]) 
vanishes. The proof of Theorem [151 relies on a non-trivial "kernel-splitting" argument that 
is discussed in Section [9] 

4. CUI AND FREEDEN KERNEL 

For § 2 , Cui and Freeden [5T] studied the kernel 

°° 1 fl 

(35) A- CF (x, y): = l + Jfl^Y) P/(x • y) = 2 - 2 log (l + J ~ * ' Y ) , x,y£§ 2 . 
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It was observed in [IDJ that this is a reproducing kernel for H 3//2 (S 2 ) as can be seen from 
the above Laplace- Fourier series expansion in terms of Legendre polynomials; cf. Section [2j 
Since the constant term in the series expansion fl35|) is 1, Proposition [TT1 asserts that the 
corresponding worst-case error for HI 3 / 2 (S 2 ) equipped with this kernel is 

/ 9 N N r— ; 

(36) wce(g[X JV ];e 3 / 2 (S 2 ))= l-^^log(l + J ' | ' 

\ j=l i=l 

The right-hand side, up to a constant factor, is known as the Cui and Freeden (CF) 
discrepancy of Xn- Note that by Theorem [121 sequences of N -point configurations that 
minimize the CF discrepancy are sequences of QMC designs for H 3//2 (§ 2 ). This fact can also 
be seen without appealing to results for spherical designs by applying the independently 
derived Theorem [H2 in the next section and using the equivalence of norms J3 In [2TJ this 
discrepancy has been used to test for uniformity of a sequence of iV-point configurations. 
Furthermore, the CF discrepancy was used in [ID] for analyzing quadrature properties of 
so-called extremal or Fekete points (these are sets of (t + l) 2 points on S 2 that maximize the 
determinant of the interpolation matrix for polynomials of degree t). Numerical data in 
[ID] suggested that for H 3 / 2 (§ 2 ), the CF discrepancy for spherical t-designs obtained using 
Fekete points as starting points, decays like 0(t~ 3 ^ 2 ); this in turn led to the discovery of 
Theorem [TJ 




5. Generalized distance kernel 



In the following we make use of the identity 

|x-y| 2 = 2-2x-y, x, y e S d . 

Reproducing kernels for H s (S d ) for s > d/2 can be constructed utilizing powers of 
distances, provided the power 2s — d is not an even integer. Indeed, it is known (cf., e.g., 
[29] ) that the signed power of the distance, with sign (— 1) L+1 with L:=L(s):=[s — d/2\, 
has the following Laplace-Fourier expansion: 



(37) (-1) L+1 |x - y\ 2s ~ d = (-l) L+1 V d ^ 2s (S d ) + «J S) Z ^ Q P f(x • y), x, y G 

t=i 

where 

V d -2 8 [b ).= / x-y da d (x)da4y) = 2 — l=r7T7o i S ' 



id 



a«:=^ (S -)(=0^. i>l. 



(3 

» _ 

(d/2 + s) 

(s) 

From these formulas one can verify that all the coefficients a\ are positive for t > L + 1 
and alternate in sign for £ < L. Furthermore, the ct^'s decay with the rate required for 



^An even more direct proof can be given by applying 1341 Theorem 2.2]. 
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coefficients in the Laplace- Fourier expansion of a reproducing kernel for HP(§ d ) as can be 
seen from the asymptotic expansion 

(39) a« ~ 2 ^-i r((d+l)/2)r( g ) 

[6J) ° e ^\{-l)^T{d/2-s)} as ^°°- 

Thus, by modifying if necessary some of the early coefficients, one can derive a reproducing 
kernel for W(S d ) for s > d/2 and 2s - d not an even integeio (cf. Section 12.41) . 

Case I. For d/2 < s < d/2 + 1 (in which case, L(s) = 0), only the constant term in f )37|) is 
negative and thus by adding any constant larger than Vd-2s{^ d ), say 2Vd_2 S (S d ), we obtain 
the following reproducing kernel for H s (S d ) which we call the "generalized distance" kernel: 

(40) K^(x,y):=2V d _ 2s (S d ) - |x-y| 2 *- d , x,y G § d . 

In particular, for s = (d+ l)/2 we get the "distance kernel" for H^ +1 ^ //2 (§ d ): 

(41) K dist (x,y):=^« d+1)/2) (x,y) = 2V_ 1 (§ d ) - |x - y| , x,y G S d , 

which for d = 2 is equivalent to the Cui and Freeden kernel in the sense that there exist 
positive constants c and C such that cifdist(x, y) < _Kcf(x, y) < C i^dist(x, y) for all 
x,yG§ 2 . 

With respect to the K^J kernel, the worst-case error for d/2 < s < d/2 + 1 is, from 
Proposition [TT], given by 

( N N \ 1 / 2 

j=i i=i j 

According to Theorem [T2l for d/2 < s < d/2 + 1, minimizing the right-hand side above, or 
equivalently, maximizing the sum of generalized distances, yields QMC designs for H s (S d ). 
This fact can also be established without appealing to any properties of spherical designs 
(and hence is independent of Theorem [8]). Indeed, Wagner [H], extending a result of 
Stolarsky [32], showed that for d/2 < s < d/2 + 1 there exists a sequence of iV-point 
configurations {x* ^ . . . ,x^ iAr } and a positive constant rj s ^ such that 

1 N N 

(43) V d ^ d ) - W2 J2H l x ,V " KM 23 '" < V.4 N ~ 2s/d > N > 2 - 

j=i i=i 

(This fact also follows immediately from Theorem |2~41 below dealing with randomized equal 
area points on S d .) Consequently, we have provided an independent proof of the following 
result. 

Theorem 16. Given s G (d/2, d/2 + 1), a sequence of N -point sets that maximize the 
generalized sum of Euclidean distances (171) is a sequence of QMC designs for H s (S d ). 



*In the case of 2s — d is an even integer, the expansion (|37p terminates after finitely many terms and so 
the Q!g s do not satisfy the appropriate asymptotic behavior (fig)) for H s (S d ). 
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We remark that for s > d/2 + 1, iV-point configurations on E> d with maximum general- 
ized sum of distances (without further restrictions) will have a limit distribution that is 
concentrated in two opposite points on E> d (Bjorck pm Remark 1 following Theorem 7]) 
and, clearly, do not lead to QMC designs. 

Case II. For s > d/2 + 1 and L as defined as above (so that d/2 + L < s < d/2 + L + 1), 
the representation ( 1371) gives rise to a reproducing kernel for IF(§ d ) of the form 

(44) /jrW(x,y):= (l - (-1) L+1 ) ^ 2s (§ d ) + Q L (x-y) + (-l) L+1 |x-y| 2s ^, x,y £ S d , 
where Ql is a polynomial of degree L > 1, 

L 

(45) Q L (x • y):=Y, (("1) L+1 "' - l) «? } Z(d, i) P £ (d) (x ■ y), x, y £ § d , 

£=1 

that simply changes the signs of the negative coefficients af \ i > 1, in (1371) . As a conse- 
quence of Theorem [121 we obtain the following. 

Theorem 17. Given s £ (d/2 + L, d/2 + L+l), where L is a positive integer, the sequences 
of N -point sets that minimize the worst-case error 

wce(Q[X N };U s (§ d )) 

1/2 



(46) 



/ TV TV 

\ i=i i=i 



Xj Xj 



i2s-d 



> i+1 VW§ d ) 



/orm a sequence of QMC designs for HP(S> d ). 

For fixed L > 1, one can avoid the introduction of the Q^-energy term in the worst-case 
error formula above by restricting attention to node sets X^,l = { x i,l> • • • > x tv,l} that are 
spherical L-designs; i.e., satisfy 

TV 

(47) Y ^J,l) = 0, l<£<L,l<k< Z(d, £). 

j=i 

For such sequences the worst-case error formula reduces to 

TV TV 

[wce(g[X^];ff(§ rf ))] 2 = ^^(-l)^i| Xj - L - Xi , i | 2 '- d - {-l) L+l V d . 

j=i i=i 

Thus, spherical L-design configurations that minimize 

TV TV 

" x ^| 2S " d ' N>c d L d 
j=i i=i 

yield sequences of QMC designs for HP(S d ) whenever 2s — d is not an even integer. 

Note that for even 2s — d = 2L = 2,4,..., the expansion (137]) terminates and is a 
polynomial of degree L in x • y. In such a case the radial (signed) generalized distance 



sl§ ) 
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(— l) L+1 |x — y| 2i log |x — y| can be used to define a reproducing kernel for H 2i (§ d ). This 
approach is explored in [IT] . 

6. Uniform distribution and low-discrepancy sequences on the sphere 

Uniform distribution. An infinite sequence (Xn) of Appoint configurations on E> d is 
asymptotically uniformly distributed on E> d if for every / G C(S d ) there holds 



For a sequence (X N ) of QMC designs for HF(§ d ), s > d/2, it follows from (125]) and 
( fj~8l) that ( 148]) is satisfied for any polynomial on S d and hence for all / e C(E> d ), since the 
polynomials are dense in C(S d ). 

Proposition 18. Given s > d/2, a sequence of QMC designs for HP (S d ) zs asymptotically 
uniformly distributed on E> d . 

Asymptotic uniformity of point sets is a relatively weak property: Even if we restrict our 
attention to / e W(§ d ) for all s, / non-constant, it is possible to construct asymptotically 
uniformly distributed sequences of node sets so that the convergence of the quadrature 
error is as slow as one likes. (Without loss of generality, we can restrict attention to 
smooth functions / whose integral over §> d is zero. One can then assign too many points 
to regions where / is positive and too few where it is negative, correcting the imbalance 
as A" — > oo as slowly as one wishes.) 

Note that asymptotic uniformity of a sequence of point sets does not imply that the point 
sets have Property R (as one can always add order o(N) arbitrary points to an Appoint 
configuration without changing the limit in (]4"8l ). nor need a sequence of point sets with 
Property R be asymptotically uniformly distributed (as one can always delete all points 
inside a fixed spherical cap without any loss to Property R). 

Low-discrepancy sequences. Unlike the situation for the unit cube, on the sphere E> d 
there is no single Koksma-Hlawka inequality (cf., e.g., [23]), as evidenced by the many 
competing Koksma-Hlawka like inequalities proposed in the literature (e.g. [2TJ, [T4]). The 
notion of worst-case error (Definition [TP]) provides a way to bound the error of numerical 
integration for sufficiently smooth functions /, see (126]) . 

The Sobolev space H^ +1 ^ 2 (S d ) with (1411) as reproducing kernel is special in the sense 
that the worst-case error wce(Q[X N ]; M^ d+1 ^ 2 (E d )) of a QMC rule with node set X N = 
{xi, . . . , xjv} has an interpretation as the spherical cap ^-discrepancy of Xn, defined by 



(48) 




as A^ — )• oo. 




N 



2 



1/2 



(49) 




(x N y. 



<7 d (C(x;0)) da d (x) sin0d0 
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where C(x; 9) denotes a spherical cap as defined in fl32l) . Indeed, Stolarsky's invariance 
principle [12] asserts that 



N N \V2 



j=l i=l / V ' a 

where V_i(S d ) is given in ( )38l) and 7^ is given in (|9]). Recalling equation ( 1421) . we recognize 
that the left-hand side above is the worst-case error for M( d+1 )/ 2 (§ d ) with kernel if dist (x-y) 
(cf. flH)))ll As a consequence of fl50|) and (H2l) for s = (d + l)/2, we obtain the following 
corollary to Theorem [ 



Corollary 19. Minimizers of the spherical cap ^-discrepancy form a sequence of QMC 
designs for U^ 1 ^ 2 (S d ) . 

A related concept is that of the spherical cap L^- discrepancy of an Appoint set on 
E> d , defined by 



(51) DC(X n ):=su V 



x N nc\ 



: C spherical cap in S a 



N 

Note that D£ 2 (X N ) < y/2D£ x (X N ). We shall establish in Section [9] the following. 

Proposition 20. Under the assumptions of Proposition [Tl\ given s > (d + l)/2, every 
N -point configuration on S d satisfies 



(52) wce(Q[X^];M s (§ d )) < c B>d D°(X. 



N) 



where c s ^ > depends on the M. s (S d )-norm, but is independent of N. 

Thus, node sets with small spherical cap discrepancy are of some interest with regard to 
numerical integration. We remark that a sequence (Xn) of ./V-point configurations on S d 
is asymptotically uniformly distributed if and only if D£ (Xn) — > as N — > 00 (see, e.g., 

ra. 

Beck [5] proved that there is a positive number c\ such that for any iV-point set on 
E> d there exists a spherical cap Cn C § d such that 



Cl 

t 53 ) N [(d+l)/2]/d < 



\Z N n I 



N 



N, 



and, by employing a probabilistic argument, that there exist C2 > and iV-point sets Z* N 
on E> d such that 



D^(Z* N )<c 2 



N [(d+l)/2]/d ■ 



This motivates the following definition. 



§ This connection was first noticed in [17] and rigorously proved in |15j . 
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Definition 21. A sequence (Z^) of Appoint configurations on E> d is said to be a low- 
discrepancy sequenc^ if there exists a positive number (3d, independent of N, such that for 
all Zn 

(54) Dr(Z N )<p d 



N [(d+l)/2]/d- 

Let (Zn) be a low-discrepancy sequence on E> d . Using fl52|) and Theorem [3j we see that 
for each s > (d + l)/2 

(55) < wce(Q[Z JV ];]EF(S« i )) < c M rf(^) < ' 



On the other hand, if wce(Q[Zjv]; EI s (S d )) < 1, we have by Lemma 1261 of Section [9] (with s 
replaced by (d + l)/2 and s' replaced by s) that, for each s < (d + l)/2, 

BAs d) , floe A0 s/(a!+1) 

(56) < wce(Q[Z N ];W (S d )) < B 3 (s,d) ^_J_ . 

From (l56i) we see that for every s G (d/2, (d+l)/2), low-discrepancy sequences on E> d have 
almost optimal order of the worst-case error for IF(§ d )), except for a power of log N. This 
leads to the following natural question: Do low- discrepancy sequences on S d form sequences 
of QMC designs for W(S d ) when s G (d/2, (d + l)/2) ? We remark that to the authors' 
knowledge no explicit constructions of low-discrepancy configurations on S d are known, in 
contrast to the situation for the unit cube. 

7. QMC DESIGNS ARE BETTER THAN AVERAGE 

As we show in this section, Theorem [7] is a consequence of a more general result dealing 
with the average value of the Appoint energy defined by a positive definite kernel on the 
sphere. Given a sequence (a^>i with a £ > and J2'iLi a £Z(d,£) convergent, we set 

oo Z(d,£) oo 

(57) A(x-y):=A{x,y) = J> Jl Y t>k {x)Y u (y) = ^ a e Z(d, t)Pf{x- y), x,y G S d . 

l=\ k=l 1=1 

Then A is a positive definite kernel on S d in the sense of Schoenberg [3?]; indeed, from 
the addition theorem, for all Appoint configurations X N = {x 1; . . . , x^} C S d and all real 
numbers a±, . . . , o>n, 



N N oo Z(d,i) 

3=1 i=l 1=1 k=l 



N 



j=i 



2 



> 0. 



"This terminology is short for low spherical cap h^a-discrepancy sequence. 

"In [2] it is proved that the spherical cap discrepancy of so-called spherical digital nets and spherical 
Fibonacci points is bounded by C./V -1 / 2 (for some explicit C), the same rate as for random points. 
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Consequently, the normalized TV-point energy of this kernel, defined by 

^ N N 

(58) A[X N ]:=— E A ^ " x *)' Xn = fr' C § d , 

j=l i=l 

satisfies *A[Xat] > 0. 

We are interested in the expected value of the energy (|58l) for N random points that are 
chosen independently and identically distributed with respect to the uniform measure on 
§ d ; namely 



EA[X N ]:= f ■■■ f A[X N }da d (^)---d 



which is also known as the spherical average of *4[Xtv] over all iV-point sets X N C S d . A 
straightforward computation yields the following result (see Section [9]). 

Theorem 22. Given the kernel A as in 057|) . 

(59) EA[X N ) = 

Remark. Theorem [22] generalizes the result [H] Theorem 6] . 

We now apply this result to the positive definite kernel K,^ associated with the repro- 
ducing kernel as given in (j2Ij) for W(§ d ) with s > d/2. By Proposition [II] 

1 N N 

KV[X N ] = ^EE^to = [wce(g[X^];ff(§ d ))] 2 , 

j=l i=l 

and hence by Theorem [22] we obtain that the expected value of the squared worst-case 
error is given by !C^(1)/N, from which Theorem [7] follows. Consequently, sequences of 
randomly chosen points on the sphere do not generate QMC designs for s > d/2. 

By the same token, Theorem [22] applied to A — as given in (|34j) yields the following. 

Proposition 23. Given s > d/2, for any positive integers t and N (not necessarily de- 
pending on each other) 



E 



N N 

■(«) 

j=l i=l 



JC ( t s \l) 
N 



Furthermore, ty s) (l) / /C (s) (l) as t -> oo. 

The implication of the first part of this result is that for a sequence of random iV-point 
configurations on E> d , the rate of convergence of the square-bracketed expression above is of 
order iV -1 regardless of the choice of t. We remark that, in contrast, for spherical t-designs 
the square-bracketed expression above vanishes. 

If, instead of choosing points randomly over the whole sphere, we stratify our approach 
by requiring that the N points be randomly chosen from N different equal area subsets of 
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§ d having small diameter, then on average we will obtain a sequence of QMC designs for 
H s (S d ) whenever s G (d/2, d/2 + 1). In the formal statement of this result we denote by 

diamA:= sup { |x — y| :x,yei} 

the diameter of the set A. 

Theorem 24. Let (T>n) be a sequence of partitions ofE> d into N equal area subsets Dj^, 
j = 1, . . . , N, such that diam-D^jv < c/N 1 ^, where c is a positive constant independent 
of j and N. Let = {x^, • • • ,~X-n,n}> where x^jv is chosen randomly from Dj^ with 
respect to uniform measure on Dj^. Then, for d/2 < s < d/2 + 1, 

(6°) Wf d < ^[{^e(Q[X N ];M^))} 2 } < 

where (3' > and (3 > depend on the M. s (S> d )-norm, but are independent of N, and (3 > 
also depends on (T>n)- 

We remark that such a sequence of partitions always exists; see 0, [121 E21 134"] . 
As the next result shows, the stratification strategy does not lead however, on average, 
to QMC designs for U s (§ d ) with s > d/2 + 1. 

Theorem 25. Let CDn) an d X N be as in Theorem\2^ Then for s > d/2 + 1, with 2s — d 
not an even integer)^ 

where (3 > depends on the M. s (S d )-norm and {T>^), but is independent of N. 

Other concepts of randomness can be also considered. Armentano, Beltran, and Shub [1] 
study point configurations on S 2 that are derived from the zeros of random polynomials. 
This will be a topic of future research. 

8. Numerical Experiments 

8.1. Point sets. Many different sequences of point sets are potential candidates to be 
QMC design sequences for IF(§ 2 ) and some s > 1. We consider the following point sets 
X N . 

• Pseudo-random points, uniformly distributed on the sphere. 

• Equal area points based on an algorithm given in [51] . 

• Fekete points which maximize the determinant for polynomial interpolation 

• Coulomb energy points, which minimize 

N N 1 

EE 

j=l i=l 



Xj 



ft A similar result for 2s — d an even integer can be obtained using an appropriate kernel as discussed 
in Section [5] 
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I Xj X j 



• Log energy points, which minimize 

N N 

EE* 

j=i 1=1 

• Generalized spiral points (cf. [SUIT]), with spherical coordinates (9j,(j)j) given by 

2? — 1 

Zj = l- — , dj = cos -1 ^), 4>j = 1.8VN0J mod 2vr, j = 1, . . . , N. 

• Distance points, which maximize 

N N 

j=i i=i 

• Spherical t-designs with N = \(t + l) 2 /2] + 1 points. 

All the point sets that are characterized by optimizing a criterion are faced with the 
difficulty of many local optima. Thus, for larger values of N, these point sets have objective 
values near, but not necessarily equal to, the global optimum. 

QMC designs for HP (§ 2 ), for every s > 1 (s not an integer), could be calculated by 
optimizing the expression ( 146]) in terms of the generalized distance |xj — Xj| 2s_2 (including 
the low order polynomial Ql or by imposing the additional constraints that the point set 
is an L-design), as discussed in Sectional We restrict attention to the point sets listed 
above, which are available from the website |45] . Here the criterion ( 14"6"j) is used to generate 
points only for the case s = 3/2 (maximizing distance sums). In all cases, except the last, 
the number of points N was taken to be a perfect square. 

8.2. Worst case error. The worst-case errors in HP (§ 2 ) for s = 3/2 are illustrated in 
Figure [TJ For all point sets, the worst case error with s = 3/2 is calculated using ( H2"j) 
and the distance kernel. Apart from the random points (which are not QMC designs, see 
Theorem [7]) all the point sets have a worst-case error for s = 3/2 that seems to decay like 
iV~ 3 / 4 , implying that they are all QMC designs for s = 3/2. 

Spherical designs with N = 0(t 2 ) are QMC designs for HP (S 2 ) and all s > 1. From 
Section O and Theorem [TBI the distance points are provably QMC designs for s = 1.5. 
The rate of decay for equal area points fits well with Theorem [2U which established that 
randomized equal area points are also QMC designs for 1 < s < 2. Moreover, from 
Theorem [251 randomized equal area points cannot do better than this. Other than the 
distance points, it has yet to be established rigorously that the non-random point sets are 
QMC designs for s = 1.5. It is rather curious that, for the non-random sequences, the 
computed worst case errors in Figure [TJ essentially lie on the same curve. 

Figure [2] plots the worst-case errors for s = 4.5 and estimates the rate of decay by finding 
a least squares fit of the form aN~P for 10 < N < 10 4 (except for the spherical designs 
which use 10 < N < 5 x 10 3 ). For all point sequences the worst case error with s = 4.5 is 
calculated using the generalized distance kernel and (1461) (for which L — 3). As expected 
from Theorem [TJ for random points the worst case error still decays like A r ~ 1 / 2 . 
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Number of points N 



Figure 1. Worst case error for IF(§ 2 ) and s = 3/2 



All the point sets, except for the spherical designs, exhibit varying rates of decay slower 
than 0(N~ s l 2 ) = 0(N~ 2 ' 25 ), indicating that their effectiveness for equal weight numerical 
integration on HP(§ 2 ) when s > 4.5 is less than optimal; cf. Theorem [TBI 

8.3. Integrating a smooth function. One expects from the error bound ff26l) that a 
putative sequence (Xn) of QMC designs for H S (S 2 ), s > 1, will play out its full strength 
when tested with a smooth test function. Some caution is needed in the choice of this 
function in order to avoid having an integrand that is accidentally too easy. Our choice is 
the Franke function for the sphere [36] defined by 

/ 0, y, z):=0.75 exp(-(9x - 2) 2 /4 - (9y - 2) 2 /4 - (9z - 2) 2 /4) 

+ 0.75 exp(-(9x + l) 2 /49 - (9y + 1)/10 - (9z + 1)/10) 

+ 0.5exp(-(9x - 7) 2 /4 - (9y - 3) 2 /4 - (9z - 5) 2 /4) 

- 0.2 exp(-(9x - 4) 2 - (9y - 7) 2 - (9z - 5) 2 ), (x, y, zf e § 2 , 
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10 1 ^™^r 




Number of points N 

Figure 2. Worst case error for HF(§ 2 ) and s = 4.5 

which is in C°°(§ 2 ) and for which 

/ /(x)da 2 (x) = 0.5328652500843890 .... 

As / G HP(§ 2 ) for all s > 1, the integration error for a particular sequence (Xn) of QMC 
designs with s* given by (jSJ), must decay at least as fast as O (iV _s * //2+e ) for any e > 0. In 
Figure [3], the faster than algebraic decay of the numerical integration error for spherical 
designs is apparent. 

8.4. Estimating s*. In order to estimate the value of s* defined by (0) we need to calculate 
the worst-case error for HP (S 2 ) with a value of s > s* (cf. Theorems 141 and [T3"|) . Similarly, 
if the error for the Franke function is decaying approximately like N~ s ^ 2 , then that point 
set can only be a QMC design for s < s. 
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Number of points N 

FIGURE 3. Numerical integration errors for the Franke function 



Table 1. Estimates of s* for d = 2 



Point set 


s* 


Fekete 


1.5 


Equal area 


2 


Coulomb energy 


2 


Log energy 


3 


Generalized spiral 


3 


Distance 


4 


Spherical designs 


oo 



Some conjectured value of s* are given in Table [fl based on the results in Figure El 
similar experiments with different values of s, and the results in Figure [31 For example, 
the equal area points have a estimated rate of decay jV~ ' 95 for the worst case error with 
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s = 4.5 in Figure [21 while the error decays like jV~ - 96 for the Franke function in Figure [3], 
leading us to conjecture that s* « 2. 

Determining the precise value of s* is very much an open question. 

9. Proofs 

Throughout this proof section we use the shortened notation Qn for a QMC rule Q[-^V| 
defined by a node set X N = {x 1; . . . , x^} C § d . 

9.1. Proofs of Section [T] results. The proof of Theorem H] requires the following lemma. 
Lemma 26. Given s > d/2, if wce(Q N ; H s (§ d )) < 1, then 

(61) wce(Q N ;M s '(S d )) < c(d,s,s') [wce(Q N ;M s ^ d ))] s ' /s , d/2 < s' < s, 

where c(d, s, s') > depends on the norms for IF(!§ d ) and H s '(S d ), but is independent of N. 

Proof. Relation (16T1) follows from Theorem 3.1 in Brandolini et al. [13]. For the sake of 
completeness we give here a proof of (161]) along the lines of the proof of Theorem 3.1 
tailored to our needs and specifically to the case of spheres § d . Throughout this proof we 

(s) 

use the canonical kernel KcJa- 

Writing 1/(1 + Xi) s in terms of a Laplace transform (cf. [H Eq. 5.9.1]), 



1 



(i + \ e y r( s 'j Jo 

and applying Proposition [TT1 we obtain for all d/2 < s' < s 

r i 2 i r°° 

(62) wce(Q N ;M s '(§ d )) = — / r s "V^t) dr. 



Here 



z(d,e) 



^(r):=^(r;x 1 ,...,x J v):=^e XeT ^ 



N 



v 2 j v v 



AT2 

j=l i=l 



1=1 k=l 

where calligraphic % denotes the heat kernel with the constant term removed: 

oo Z(d,Z) oo 

tf(T,x,y):=5>-^ £ ^, fc (x)y <jfc (y) = £ e"^ P, (d) (x • y), x,y G § d , 

£=0 fc=l £=0 

which is the fundamental solution to the heat equation du/dr + A* d u = on M + x S d . 
Interchanging integration and summation in (I6"2"l) is justified, because the heat kernel is 
uniformly continuous on [e, oo) x § d x S d for e > 0. 

Let d/2 < s' < s and set e: = [wce (QAr;ff(§ d ))] 2 / s . (Then e < 1 by assumption.) We 
split the right-hand side of (162"1) into three parts and use the following estimates: (i) For 
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"large values" of r, 
1 



IV) Ji 



OO I 
S'-I-T, ' - ' 



r s e T g{r) dr < 



r^) r( s ) J 



r s - l e- T g{r)dr 



IV) r(s') 



where we used fl62l) with s' replaced by s. 
(ii) For "medium values" of r, 



< (e/2)' 



T -V'-V' S (r)d r <<^ 

e/2 1 I s J 

-.r(a) i 



e/2 



r(*0 r( s ) 7 



00 r(?) 



r( s ') 



r( S ' 



-e . 



(iii) For "small values" of r, we appeal to the small time Gaussian estimate on the heat 
kernel (cf. [E]); that is, if < r < e/2, then for some c > 

T d / 2 H(T,x,y)<ce d / 2 H(e^,y), 

which in turn implies, on the assumption that g{r) is uniformly bounded on [0, 1), 



N N 

r d/2 9 (r) = ^ E E rd/2H ^ *0 - r " /2 ^ C£d/2 (9(e) + 1) < c'"e d/2 . 
i=i i=i 

Postponing the proof of the uniform boundedness, it then follows that 



r< 



< 



r(«o 



< 



r d/2 



r d/2 



e/2 



. s '- d /2-l e -r dr 



IV) 



b/2 



r «'-d/2-l d 



r 



r(s') s'-d/2' 



From (i), (ii) and (iii) we get the required estimate 



wce(Q7v;H s '(§ d )) 



< 



r(«) + 2 s- s > m , c" 



r(s ! ) r(s') r(s')s'-rf/2 

< c"" [wce(Qjv;H s (S d ))] 2s ' /s . 
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It remains to prove uniform boundedness of the g(r) for < r < 1. (We use (|26|) and 
the fact that J §d 7i(r, x, x fe ) dcrrf(x) = 0.) 

1 N 1 N 1 N 

< 9(r) = ^J2n T, u ^ x ^> x «) ^ wce (^; ff ( §d )) ^ I>( r > ■» x ^) 

3=1 1=1 



wce(Q N ;U s (§ d )) 



wce(Q N ;U s ($ d )) 



< wce(Q N ;M s (S d )) 



oo Z(d,£) / N 



ifc (xj) >W 



£=1 k=l 



z{d,e) 



i=i 



2A«t 



=1 k=l 



1 * 



1 1/2 



£=1 fc=l 



1 " 



1/2 



sup{(l + A,) s e- A ^} 



[wce(g iV ;e s (§ d ))] 2 sup{(l + A £ ) s e-^ r } = r s sup {(1 + \ e ) s e~^ T } . 



£>1 



£>1 



The function (1 + X) s e Ar has a unique maximum atA = s/r— 1 with value s s r s e T s . 
Hence 



< g(r) < s s e T - s < s s e 



s _1— s 



□ 



Proof of Theorem^ If (X^) is a sequence of iV-point QMC designs for IfF(§ d ) for s > d/2, 
then wce(Q7v;H s (§ d )) ^ as iV ^ oo. Thus wce(Q N ; H s (§ d )) < 1 for all N > N for 
some iVo > 0. By (|6T|) 



wce(Q 7V ;H s '(§ d )) < c(d, s, s') [wce(Qjv; H s (§ d ))] < c(d,s,s' 



c(s, d) 
N s l d 



s'/s 



]S[s'/d 



for every N > N for each d/2 < s' < s. The finitely many exceptions with N < N 
satisfy the last inequality with a possibly larger constant c", depending only on the norms 
of M s {S d ) and U s '(§ d ). Consequently, (X N ) is a sequence of QMC designs for U s ' (E d ) for 
d/2 <s' <s. □ 

Proof of Theorem® Let d > 2. By Theorem [8] there exists a sequence (Y Nt ) of spherical 
t-designs Y/v t with N t :=m ( it d points (t > 1) for some suitably large positive integer 
Furthermore, the theorem states that there exist spherical t-designs for every cardinality 
> rrid t d . Thus we can fill the gaps in the sequence (Vjv t ) by adding spherical t-designs with 
N points for N t < N < N t+ i. If necessary we choose for N = 2, . . . , spherical 1-designs 
with N points; that is, configurations with centroid 0. This gives a new sequence (Yw)./v>2- 
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By Theorem [1] there exists a constant C(s,d) > such that 



(63) sup 

feH s (S d ), 
II/IIh»<i 



N 



„ C{s,d) _ C(s,d) (N^X 

- j~ s ~ N s/d y i J 



for all N > 1. Since 



Crf = —7— < — ^— < — ^— = c d ' — ^— for all iV t < iV < N t+1 , 
the right-hand side of fl6"B"j) satisfies for all JV > 1 

« < fejHl±W < ^M, c - (M):=2 vf c(M). 

]V s / d \ t J ~ N s / d ~ N s l d 

Consequently, (Yjv)jv>2 is a sequence of generic QMC designs. □ 
9.2. Proofs of Section [3] results. 

Proof of TheoremUM For N > 2 let = {x* N , . . . ,x* N N } be as in Theorem [TJ] and 
Y/v = {yi,N, • • • 7 Yn,n} be as in Theorem [6j By minimality of the X^s, for every N > 2 

E E k w (4* • X U) ^ E E • y^). 

J=l 8=1 J = l 8=1 

Hence, by Proposition [Til 

wce(Q[X^];ff < wce(gp^]; < 
where the last inequality follows by Theorem [61 □ 

For an iV-point configuration on S d the worst-case error formula (T29]) can be ex- 
pressed in terms of the polynomial part (I34p and the tail part 



(64) 4 ) Z(d,£)Pi d) (z), z e [-1, 1]; 

e=t+i 

that is, 

. JV n n 

(65) [Ke(Q[^]; W (S d ))} 2 = — £ E • x *) + N2 E E tf'fo ' x *)- 

j=l i=l j=l i=l 

Since and are positive definite in the sense of Schoenberg, both double sums 
in (1651) are non- negative (cf. Lemma [9]). Consequently, a sequence of QMC designs must 
satisfy the following necessary conditions. 
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Proposition 27. For a given sequence (X N ) of QMC designs X N = {x 1)A r, . . . ,it NtN } for 
M. s (S d ), s > d/2, there exists c(s,d) > independent of N such that, for all integers t 

1 N N XI JM2 ! JV JV r / a 12 

2^ 2^ ^ x ^ ' Xi .^ - N 2s/d ' 2^ 2^ x * ^ x ^ ' Xi >^ - ]v 2s / d ' 

j=l i=l j=l i=l 

As discussed in [TjSj, one cannot expect to get the desired estimate 0(t~ 2s ), which trans- 
lates into (D(N~ 2s / d ) if N x simply by term wise estimation of the energy corresponding 
to the tail portion < xf > of the reproducing kernel. The papers [161 [23 I2H] thus devised a 
method of splitting the reproducing kernel into a polynomial part, which could be inte- 
grated exactly by an appropriate numerical integration rule, and a remainder part small 
enough to be bounded in a simple way. In [17] this method is applied to derive the following 
estimate for the worst-case error for a sequence (X^) that satisfies Property R. 

Theorem 28. Under the assumptions of Proposition [77], if (X^) is a sequence of N- 
point configurations on E> d that satisfies Property R, then there exists a constant c > 
independent of N such that for integers t with t x N l l d 

N N 

[wce(Q[X N }; 2 < A E K^C** • x *) + 0{N" 2s / d ) as N oo. 

i=i i=i 

Proof of Theorem[T5[ Let (Xjy) be a sequence of Appoint configurations on S d satisfying 
Property R. If (Xjy) is a sequence of QMC designs for IF(§ d ), then by Proposition E7J for 
all integers t 

N N 

j- 2 E E ^(4 • x ') = 0(N~^ d ) as AT oo. 

In the converse direction, if (Xjy) satisfies the above condition, then by Theorem [28] 

[wce(Qjv; ff(§ 2 )] 2 = 0{N~ 2s/d ) as A^ oo, 
and hence (Xjv) is a sequence of QMC designs for HP(§ d ). □ 

9.3. Proofs of Section [6] results. 

Proof of Proposition IWl We recall that two different kernels for the same EI s (S d ) yield 
worst case errors that can differ by at most constant factors, and also that for the particular 

(s) 

case of the canonical kernel -fQan, the worst case error decreases monotonically with s for 
s > d/2, cf. relation (}3~Tj) . Using these facts and (150]) . which gives the worst-case error 
with respect to -Kdist, we have for s > (d + l)/2 the following estimates: 

wce(Q[X N ];W($ d )) < (3 l wce(Q[X N ]; M^W 2 (S, d )) < 2 D^(X N ) < (3 3 D°JX N ), 

where the positive constants 0i, 02 and 02 depend on the chosen norms. □ 
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9.4. Proofs of Section \7\ results. 

Proof of Theorem\2^ Let A L (z) be the truncated series 

L 

A L (z):=J2cLiZ(dJ)Pi d \z), ze [-1,1]. 
i=i 

On separating the diagonal and off-diagonal terms of the double sum 

N N 

x _ r„ „ 1 r- 



A L [X N ] 



N 2 



{Xi, . . . ,Xat} C 



j=l i=l 



we obtain 



EA L [X_ 



N 



N 



N N 



i=i i=i 



d(T d (xi) • • -do-rffxArl 



+ mEE / / ^ L ( x i • x *) d ^( x i) d ^( x * 

7V ~[~iJndJ^ 



A L {1) N(N-l) 



N 
N 



+ 



iV 2 



L Z(d,£) 
£=1 k=l 



3^ fc (x)dcr <l (x] 



For constant coefficients 1 this is the result of Theorem 6 in [4T1. 

The quantities ^[Xjy], ^4. 2 [Xjv], • • • form a point-wise non- decreasing sequence of non- 
negative (o~d ■■■ a,i) -measurable functions with limit function ^4[Xjy]- By the monotone 
convergence theorem it follows that 

EA[X N ] = [ ■ ■ ■ I A[X N ] d<r d (xi) • • • da d (x N ) 

A L (1) A(l) 



lim EA L \X N ] = lim 



This completes the proof. 



□ 



Proof of Theorem^24\ We follow the proof idea leading to [3U Theorem 2.2]. Let T>^ = 
{Dj } N, . . . , Dn,n} be an equal area partition of E> d into subsets with small diameter; that 
is: Uf^Dj^ = E> d , where o~d(Dj^ H -Dfe,v) = for all j, fc = 1, . . . , N with j ^ k and 
cr,i(Dj t N) = 1/N; furthermore, diam-D^jv < c/N l / d for some c not depending on N. Each 
Dj,N is equipped with the probability measure 



(66) 



0~d 



fJ>j,N- 



D 



j,N 
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Let d/2 < s < d/2 + 1. Then the expected value of the squared worst-case error 
[wce(Q[Xjv];H s (§ d ))] 2 for the space W{E> d ) provided with the kernel K$ given in gOJ) 
when the j-th node is chosen randomly from -D^jv (with respect to uniform measure on 
Dj jN ) is given by (see ( 1421 ) 



E 



{wce(Q[X N ];M s (S d ))f 

= ! ■■■ I [wce(g[X^];e s (§ d ))] 2 d / i 1 ^(x 1 )---d/i^(x^) 



N N 



v d ^ d ) 



N 2 EE 



I X j 



1 2s— d 



d/i J - A r(x i ) dfi i)N (Xi 



j = l i=l ° D 3,N J Di,N 



/ / |x-y| 2s d da d {x) da d (y) 
J§ d Js d 

1 ^ 

/ l x -yl 2s_dd /W x ) d /^vv(y) 



Since the first double integral in brackets equals V^„ 2 s(S d ) (see 

N 



we deduce that 



E 



{wce(Q[X N );W($ d ))} 2 =^Y1 



N 2 



x-y| 2s d d/i jiA r(x) dfi j;N (y) 



j = l J D J,N JDj >N 

N 



N N 



i=i j=i 
The lower bound in (160 p follows from Theorem [31 



□ 



Proof of Theorem{25[ Let rf/2 + L<s<<i/2 + L + l for an integer L > 1. Arguing as in 
the proof of Theorem l24"t but using the kernel in (I44p . we obtain 



E 



{wce(Q[X^] ; e s (§ d ))}' 



N 2 



N 

E, , 

j=l J Dj, N J Dj M 
N 



[Ql(1) - Ql(x ■ y)] d/i jVV (x) d/ij )JV (y) 



/ / (-l) L+1 |x-y| 2s - d d^(x)d^(y) 

iV j=i J Dj,n JD J>N 



where we used the fact that the Laplace-Fourier expansion ff45]) only contains Legendre 
polynomials with £ > 1 and thus 



N N 

I I QL(x-y)d/i^(x)d^(y)= / / Q L (x-y)d^(x)da d (y) = 0. 

!V ;=i i= i "'•Di.JV ^ A,iV ^ ^ 
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For L = 1 the definition of Ql given in ( )45|) and the fact that P\ (x) = x yields 

Qi(l) - Qi(x • y) = -a[ s) Z(d, 1) (2 - 2x ■ y) = -a[ s) (d + 1) |x - y| 2 , 

where < by (I38p . More generally, using the following hypergeometric function 
relation for the polynomials P^ d \x) (see, e.g, [TJ Eq. 18.5.9]), 

we have the following representation for Ql(x • y) in terms of even powers of distances: 
Qx(l)-Qx(x-y) 

- &(D - E E «-D i+I -' - 1) ') ^^ (^F)' 

p =o ^ ' 'p^' ^ ' 

= -EjE((-^-i)^(*<)(^ + --i),} w ( 



i x -y| 



For small distances the dominant term is the square of the distance. Since ( — 1) L+1 ^af 
is positive by (I38p . it follows that the coefficient of (|x — y|/2) 2 , 

is positive and therefore 

Ql(1) - Ql(x ■ y) = fi[ s) P^p^j 2 + 0(\x - y| 4 ) as |x - y| 0. 



Hence, 



E 



{wce(Q[X^];ff(§ d ))} 2 



(67) > i JL 



^2 / |x-y| 2 d/i iiA r(x)d/i ii 7v(y) -"^tv, 



4 iV 2 

3 

where (as N — > oo) 

JV JV 

7^ = C(= £ [diamD^] 4 ) + 0(— £ [diamD^; 
l by J i=i i=i 

= OiN-V*- 1 ) + 0(N~ 2s/d ). 



2s-d\ 
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Next, we observe that for any c" > the following inequalities hold: 



/ / |x-y| 2 d/x ii j V (x)d// i) jv(y) > / / |x - y| 2 dfi jiN (x) dn j)N (y) 

\x-y\>c"/N 1 / d 

>(c") 2 N-V d ! ! d^ A (x)d^ A (y) 
(69) Jd hn Jd,. n 



\x-y\>c" /N 1 /* 



z 'f N -2/dL_ f f d ^ (x)d ^ (y ) 

JD 4 N J D-i AT 



\x-y\<c"/N 1 / d 



Since fj,^ is a probability measure on -D^tv and ad(Dj^) = l/N, we can bound the above 
double integral by 



/ / d/i jVV (x)d/i iiAr (y) < / / d/i jiA r(y)d/i iiAr (x) 

JD,. N JDi. N J D jtN JC(x;6') 

H jjN (D j:N ) fi jiN (C(x; 9')) 



(70) Ix-y^c'ViVVrf 

a d {D j)N ) 

where 2sin(6''/2) = c" /N 1 ^. An application of the Funk-Hecke formula gives (cf., e.g. 

TO 

a d (C(x; #)) = [2 sin(^/2)] d {l + 0([2 sin(fl'/2)] 2 ) 

Hence, for N sufficiently large, 

ad(C(x; 6 ' ] =N l -^l U*/ N V*) d U + 0((c"/N 1 / d ) 2 )) < 2-^ (c") d . 
Od{D j:N ) d u d I J d uj d 

By fixing c" (which now depends only on d) to be sufficiently small, we can always achieve 
that the double integral on the left-hand side of (1701) is bounded from above by 1/2. 
Therefore, for sufficiently large N and j = 1, . . . , N, we deduce from (R)9"|) that 

/ / |x- y | 2 d/i^(x)d/i^( y ) > ^N~ 2 ' d . 

JD j>N JD jiN z 

Combining this estimate with ( 1ST)) and (168)) we obtain for s > d/2 + 1 and 2s — d not an 
even positive integer 



E 



{wce(g[X JV ] ; e s (§ d ))} 2 ] > ^ + oiN- 4 ^ 1 ) + o(N- 2s ' d ) 

> f3N^ 2/d ~ 1 = f3 N~ 2W2+1)/d , 



where the positive constant /3 depends on the H s (§ d )-norm and the partition sequence 
(T> N ), but is independent of N. □ 
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